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Abstract 

o . 

, In this paper we introduce a PDE system which aims at describing the dynamics of a dispersed phase of 

particles moving into an incompressible perfect fluid, in two space dimensions. The system couples a Vlasov-type 
equation and an Euler-type equation: the fluid acts on the dispersed phase through a gyroscopic force whereas 

Qthe latter contributes to the vorticity of the former. 
' First we give a Dobrushin type derivation of the system as a mean-field limit of a PDE system which describes 

i the dynamics of a finite number of massive pointwise particles moving into an incompressible perfect fluid. This 

last system is itself inferred from the paper [17], where the system for one massive pointwise particle was derived 
as the limit of the motion of a solid body when the body shrinks to a point with fixed mass and circulation. 
Then we deal with the well-posedness issues including the existence of weak solutions. Next we exhibit the 
• Hamiltonian structure of the system and finally, we study the behavior of the system in the limit where the 

' mass of the particles vanishes. 

?2 ' 1 Introduction 



> 



1.1 The model at stake 

This paper is devoted to the following PDE system 



dfUJ + divxi^u) = 0, (1) 

at/ + div,(/e) + 

PsJ ■ where 



9t/ + div,(/e) + div5(/(e-u)^) = o, (2) 



u ■= K[uj + p] and p -.^ [ fd£,. (3) 

. ^ , This system describes a dispersed phases of particles (or spray) moving into an incompressible perfect fluid, using 
' a kinetic/fluid approach. The fluid is described via macroscopic quantities : uj stands for its vorticity and u for 
^ [ its velocity, both depending on the time t G IR+ and the space variable x S M^. The spray is represented by a 
■ - - I probability density function : /, depending on t, x and ^ € R^, the velocity of the particles. The notation K in 
^ stands for the Biot-Savart operator associated to the full plane; it maps a reasonable scalar function g to the 
vector field 

K[g]{x) := / H{x-y)g{y)dy, (4) 

where 

H{^) = (5) 
2tt \x\- 
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Here the notation x'^ stands for a;"*" = {—X2,xi), when x = (xi,X2). Let us also observe that 

H{x) =W^G{x). 

where G is the 2d Newtonian potential 

G{x) ■■=^H^) 

and V-^ denotes the vector field {—82, di). Moreover K[g] is divergence-free and such that curli^[5] = g. Eq.© 
defines the vector field u which describes the fluid velocity. This vector field u couples the Vlasov equation ([T]) and 
the Euler equation ([2]): the fluid gives a lift acceleration to the particles (the divj term in ([2|)), whereas the spray, 
as a whole, contributes to the vorticity of the former (equations Following the naming of (28], let us stress 

that the system (H])-® describes the behavior of a thin spray: there is no interaction (collisions, coalescences) 
between particles thus the kinetic equation is linear in / (given u) and furthermore the volume fraction occupied 
by the spray is neglected. 



1.2 A comparison with some other fluid/kinetic coupling 

Numerous examples of fluid/kinetic coupling have been studied since the seminal models introduced in jlOl 130 ) . 
but the coupling is usually made via a drag force and not a lift force. Let us compare our system with some of 
the recents models introduced in the literature. An alternative way to describe the system H])-© is to use the 
following velocity formulation: 

dtu + divx{u ® u) + Vp = {pu~i)^, (6) 
divj; u = (7) 
9t/ + div,(/0+div^.(/(^-K)^) = 0, (8) 

where 

p := [ fd^ and j := / f^d^ (9) 



Formally one can obtain the previous formulation from the system ([J)- ([3]) as follows: one integrates the equation 
([2]) with respect to ^ G to obtain the macroscopic conservation law: 

dtp + div, J ^0. (10) 

Then one adds the equation pO|) to the equation ^ to get 

dt{uj + p) + divxij + uju) ^ 0. (11) 

But, taking into account that u is divergence free, 

dt{uj + p) + divj, ((cj + p)uj — curl (^dtu + divx{u (8) w)) and diva;(j — pu) — — curl [{pu — j)^), 

so that (fTTjl now reads: 

curl (i9tu + div2;(?/ ® u) — {pu — j)^) = 0, 

what means that dfU + div^iu ^ u) ~ {pu ~ j)^ is a gradient, thus yielding 

The system (7) — (9) of [5] reads, with the present notations and for a homogeneous incompressible fluid, 

(?tu + div2;(u ® u) + Vp = j — pu, (12) 

div.u = 0, (13) 

9J + div,(/0 + div5(/(u-C)) = 0, (14) 

This can also be seen as the model (1) — (3) of [TT] where the diffusion term of the kinetic equation has been 
dropped. 

The model p2|) - p4)l involves the drag force (j — pu) whereas the model ([6])- ([8]) involves the gyroscopic force 
— {j — pu)^ . Let us mention that if several papers have dealt with spray models with drag effect, see for instance 
the introduction of |llj for more references, it is, to our knowledge, the first paper where a spray with a gyroscopic 
effect is considered. 
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1.3 Plan of the paper 

In the next section we recall the basics from optimal transportation theory that we will need in the sequel, in 
particular in Section [3] where we will derive the system (HJ-Q as a mean-field limit of a PDE system which 
describes the dynamics of a finite number of massive pointwise particles moving into an incompressible perfect 
fluid. This last model originates from the paper [T7] and can be seen as an extension of the point vortices-wave 
system introduced by Marchioro and Pulvircnti to the case where the vortices have a non vanishing mass. We will 
obtain ((J)-© in the mean-field limit in a regularized setting by adapting the famous approach used by Dobrushin 
[15] for a regularized Vlasov-Poisson system. 

In Section m we will give some existence results for the solutions of ([I]) -([3]) which include the case of quite 
irregular initial data. 

In Section [5] we will give a uniqueness result for the solutions of (II])-([3]) with w and p bounded. We will follow 
the analysis of the Vlasov-Poisson system performed by Loeper in [53] which makes use of optimal transportation. 

In section [6] we exhibit a Hamiltonian structure of the equations by introducing a Poisson structure on the 
manifold V of the pairs (w, /) such that equations H])-© yields some Hamiltonian ordinary differential equations 
for any smooth functional of a solution (w, /). 

Finally in Section [7] wc study what happens when the mass of the particles immersed into the fluid vanishes. 
We will see that the particles are then convected by the fluid velocity so that the system degenerates into the usual 
incompressible Eulcr equations. We will follow a strategy used by Brenicr to deal with the gyrokinctic limit of the 
Vlasov-Poisson equations. 

1.4 A few general notations 

In the whole paper, N will denote the set of nonnegative integers including 0, whereas the notation N* will stand 
for N* := N \ 0. We will simply denote | • | the Euclidean norm in R^, where d denotes a positive integer. 

A mapping / from W'-^ to W^'^ is said Lipschitz if there exists a constant c > such that for any x,y in M'', 
|/(x) — /(y)| ^ c\x~y\. Then the smallest admissible constant c defined a norm denoted ||/||Lip and the associated 
Banach space is noted Lip(M''i , R'^^), or simply Lip(R''i) if 3,2 = 1. 

When no ambiguity arises on the time interval [0, T] considered, we will denote by L^{L^) and Lf{L'^^ ^) the spaces 
LJ'([0,T];L«(R2)) and Lp([0, T]; L«(R2 x R^)), for all p,q e [l,oo]. Similar notation for '^t°iE) for continuous 
functions of t having value in some functional space E. When no variable is added in index, lp denotes the 
norm of h in all its variables, same thing for Sobolev norms. Thus for example, if h is defined for (i, x) € [0, T] x R^, 
||/i|jLp denotes its LP{[0,T] x R^) norm, while ||ft.(t)||LP denotes the L^'(R^) norm of h{t). Given an open subset il 
of R'', ^(f2) denotes the space of test functions having compact support in ft, S)'{H) is hence the associated dual 
space of distributions. 

We will use the same notation C for some constants which may change from line to line. We will sometimes use 
Ca to precise that the constant yet depends on a. 

2 A few tools from optimal transportation theory 

We fix from now on a few notations from optimal transportation for which we refer to [U 129] . Here d (it may be 
indexed), denotes a positive integer. It will be taken equal to 2 or 4 in the sequel, as the vorticity uj{t) depends on 
a; e R^ while the density f{t) depends on {x,S) in R^ x R^. Note also that / being a probability density function, 
it is always nonnegative, whereas we do not specify any sign condition for the fiuid vorticity cj. 

2.1 Measures 

We will denote by 7W(R'^) the set of signed measures and by A^+(R'') the set of finite measures on R''. We endow 
A^(R'') with the standard narrow convergence, given by the duality with continuous and bounded functions, that 
is to say that a sequence {vn)n in A^(R'^) narrowly converges to v in A^(R'^) if for any cf) in the space '^^{K^) of 
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the continuous and bounded functions over Mf^ (taking real values) , 

4)dvn — > I (pdv. 



For a measure v in JV[{W^) we denote by \v\ its total variation measure, and by and v~ its positive and negative 
parts, all three given by the Jordan-Hahn decomposition so that v = — u~ , and \v\ = + v~ . 

We denote V{M.'^) the set of the probability measures that is the subset of M.'^{W^) verifying i'(R'^) = 1. 

Given a measure v in Al(R''i) and a map r : we define the pushforward measure t^v ^ A^(R''^) of 

1/ by r by 

t#v[B]:^v[t-\B)], (15) 

for any Borel subset B of W^'^ . 

Given a measure 7 over the product space R'* x = M?'^ and the projections tt^ and tt^ on it factors, we define 
the first and second marginals of 7 as the measures 7r^7 and 7r|t7. This means that for any Borel set B of R"^, 

T^#l[B] = i[B X R^] and 7r|7[S] = 7[R'' x B]. 
For any p £ [1, oo[ we introduce the subspace of signed measures having a finite p-th moment 

7Wp(R^) := (i^ e X(R'') : / \xY'd\v\<oo 



and defined similarly M + {W^) and VpiW^). 

2.2 Wasserstein distances for measures 

We will make use of the Wasserstein distances. Let us first recall the standard definition for measures. 

Definition 1. Let p he a real satisfying p ^ 1, i/i and V2 be in M.^{W^) two finite measures on such as 
^i(R'^) = V2{W^). We define the Wasserstein distance Wp{vi^V2) of order p between this two measures by 

Wp{vi,i'2) := (m [ -f{x,y)\x-yfdxdy] (16) 

where the infimum is taken over all 7 G A^^(R'' x R'') with marginals vi and V2. 

It can be shown that for real p satisfying p ^ \, for any positive real number a > 0, Wp satisfies all the axioms 
of a metric on the subset of M.p{W^) of measures verifying v{W^) = a. Furthermore, convergence with respect to 
Wp is equivalent to the narrow convergence of measures plus convergence of pih moments. 

2.3 Wasserstein distances for signed measures 

For signed measures, there are several ways to define some Wasserstein distances. Let us refer here to [3]. In the 
sequel we will only need to consider some cases where the positive parts have the same total mass, and the negative 
parts as well, this motivates the following definition 

Definition 2. We will say that two signed measures vi and V2 on are compatible whenever i/j^(M'') = j/^(R''). 
More generally, we will say that a family of signed measures is compatible if any pair of its elements is compatible. 

Yet we will distinguish two cases depending on whether p is equal to 1 or not. Let us start with the second 
case. 

Definition 3. Let p e]l,oo[. Given two compatible elements vi, V2 of M.p{W^), we define the Wasserstein distance 
Wp{vi,V2) of order p between these two measures by 

Wp{v^,V2) := {Wp{vt,vtf + Wp{v^,v:, YY. (17) 
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Of course we recover Definition [T] in the case of positive measures. We will use Definition [3] with p = 2 in Section 

El 

In the case where p = 1, we extend Definition [T] by setting 

Definition 4. Given two compatible elements vi, V2 of MiiW^), we define the Wasserstein distance Wi{vi^V2) 
between these two measures by 

Wi{ui,U2) Wi{u+ + v^.u^ + ut)- (18) 
Such a definition allows us to keep the well-known Kantorovitch duality : 
Proposition 1. Let vi and V2 be two compatible elements of MiiW^). Then 

Wi{vi,V2) / (t>d{vi-V2), (19) 

where the supremum is taken over the unit ball o/Lip(R'*). 

Proof. This formula is well-known for positive measures, see for instance [2l[29]. Moreover from Definition |4l we 
infer that 

Wi{v^,V2) = sup / 4)d{{v^ +v~)- {v~ +v^)), 
if) Jw 



= sup / (j}d{l^i - V2). 

□ 

Let us give a few properties of the distance Wi. The first one is that it is really a distance: this comes easily 
from Proposition [T] For instance the positivity and the symmetry are quite obvious whereas the definiteness comes 
from Riesz theorem and the triangle inequality can be checked as follows: given three compatible signed measures 
vi, V2, 1^3 in A^(R'^), one has (the supremum is always taken over the unit ball of Lip(R'')) 



Wi{i'i,v:i) = sup / (l)d{ui- Us), 

= sup / 4>d{vi~V2 + V2-Vz), 



< sup / 4>d{vi — V2) + sup / 4>dv2 - vz), 

= Wi{l^i,V2) + Wi{v2, V-i)- 

Note also that the distance Wi satisfies the following. 

Lemma 1. If vi,V2 and v^, 1^4 are two pairs of compatible elements of Mi(R'^), then the two measures vi -l-v^ and 
V2 + are compatible and we have 

Wi{vi + 1^3, V2 + Vi) < Wi{vi,V2) + Wi{yz,^i), 

with equality in the particular case where = v^, what amounts to saying that the distance Wi is translation 
invariant. 

Proof. The compatibility is straightforward. Thanks to Proposition [T] we have (supremum is taken over the unit 
ball of Lip(]R'')) 

Wi{vi+ 1^3,1^2 + 1^4:) = sup / (/)d(z^i + 1^3 - 1^2 - i^4), 

= sup ( / (j) d{i^i - 1^2) + / (j) d{i^3 - 1^4)) 
^ sup / (^d(i^i — 1^2) + sup / (\)d{vi — V4), 
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The case where I'a = z/4 is straightforward. □ 

Remark 1. Lemma]^ can also be seen as a consequence of the fact that the distance Wi can he associated to a 
norm, see UBj- 

Remark 2. Applying Lemma[J\ with {v^ , 1^2 1 ^2 ^ ^i) instead of {i^i, 1^2, fs, V4) we have that Wi{vi, V2) is less than 
Wi{i'i ,1^2) + Wi{v^ ,V2), which is the result given by the formula in Definition\^for p ~ 1. 

Lemma 2. Consider two compatible elements vi, V2 of M^\{W^^"). Let t : R'^^ — > R"*^ be a Lipschitz map. Then 
Proof. Thanks to (|19p we have 

Wi{t^Vi,T^V2) = sup / (j)d{T^Vl ~ T^V2), 

~ sup / (f) o T d{vi — V2) ^ 

f (f) O T 

= Iklkip sup / — d{vi-V2), 

<p JK<*2 ||T||Lip 

sC ||r||Lip sup / i' d{iyi - 1^2) , 

■ip JR<1 
= IItIIlIp Wl{vi,V2), 

where the index </> denotes a supremum over the unit ball of Lip(R''^) and a supremum over the unit ball of 
Lip(R''i). □ 

In the particular case where the map t of the previous lemma is one of the projection maps defined above, we 
get the following result 

Lemma 3. Consider two compatible elements vi, V2 of M.\{W^ x R''). We have, for i = 1,2, 

Wi{'k)^VI,'k''^V2) SS Wi{vi,V2). 

2.4 Wasserstein distances for vector measures 

As we will deal with the coupled system (IT])-(l2|) for which there are two unknowns: u and /, we will use the 
following extension of Definition [31 

Definition 5. Let di and d2 he two positive integers. Consider two pairs vi, U2 and ai, 1T2 of compatible measures 
(respectively in A^i(R'*i) and Mi{R'^^)). We introduce the couples := (1/;,^;) S Mi{R'^^) x Mi{R'^^), i = 1,2 
and define the associated Wasserstein distance Wi{iii, (12) between fii and fi2 by 

Wiiiii,ii2) := Wi{iyi,V2) + Wi{ai,a2). 

3 Derivation of the system as a mean-field model 

3.1 Massive vortex- wave system 

In |17] we derive the following system for the motion of one massive vortex of mass m > and circulation 7 G R 
in a perfect incompressible flow. Recalling the notation 
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the system reads 



dtuj + div4uju) = 0, (20) 
u{t,x) := K[Lo]{t,x) + -iH{x ~ h{t)) (21) 

mh"{t) ^ -i{h' {t) - K[uj\{t,h{t))Y , (22) 

w|t=o - c^o, h{{)) = /lo, h'{Q) = hi, (23) 

(24) 



This system was obtained by considering the motion of a sohd body in a two dimensional incompressible perfect 
fluid, when the body shrinks to a pointwisc particle at the position h{t) with a fixed mass m > and a fixed 
circulation 7 G M around the body. Equation (PH)) describes the evolution of the vorticity uj of the fluid: it is 
transported by the velocity u obtained by the usual Biot-Savart law in the plane, but from a vorticity which is the 
sum of the fluid vorticity and of a point vortex placed at h(t) with a strength equal to the circulation 7. Observe 
for instance that the velocity u is divergence free and can be written as u = K[w + 7^/1(4)]. 

Equation means that the shrunk body is accelerated by a force similar to the Kutta-Joukowski lift of 
the irrotational theory: the shrunk body experiments a lift which is proportional to the circulation 7 and to the 
difference between the solid velocity and the virtual fluid velocity obtained by the Biot-Savart law in the plane 
from the fluid vorticity, up to a rotation of a 7r/2 angle. Let us refer here to the textbooks of Childress [T3] or 
Marchioro and Pulvirenti [25| for a discussion of the Kutta-Joukowski force. See also Grotta-Ragazzo, Koiller and 
Oliva |18| , where they consider a similar system of a point mass embedded in an irrotational fluid and driven by 
the Kutta-Joukowski force. 

The system ([201)- (EH) can also be seen as a variant of the vortex- wave system introduced by Marchioro and 
Pulvirenti cf. for instance Th. 6.2 of [53] and [5^, and recently studied in [3T] and in [7]. Actually the system 
reduces to the vortex- wave system when the mass m is set to 0. 

3.2 particles 

Let us now generalize the previous system to the case of N pointwisc particles of mass to^, of circulation 7^ and of 
position hi{t), for i = 1, •..,-/V, moving into a perfect and incompressible planar fluid: 

dtoj + divj. (wu) = 0, (25) 

N 

u{t,x) = if M(t, x) + Y,ljH(x - h,{t)), (26) 

mX{t)=j,(h[{t)-v,{t,h,{t))y, (27) 
v,{t,x) = K[uj]{t,x) + J2ljH{x - /ij(t)), (28) 

w|t=o = wo, h,iO) = /i.,o, KiO) = h,,i. (29) 

Let us observe that in (j28p the self-interaction is omitted since the index in the sum runs only over j ^ i- The 
derivation of the system ((25|) -([29 |) from the motion of A'' solid bodies in a two dimensional incompressible perfect 
fluid, when the bodies shrink to pointwisc particles, is the object of a paper in preparation. 

Once again if one sets all the masses equal to one recovers the vortex-wave system of Marchioro and 
Pulvirenti. 



3.3 Mean-field limit 

We want to study the mean-field limit of the system (p5)) - (p9)) . that is the limit system obtained by the empirical 
measure 



1 ^ 
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when N goes to infinity, with an appropriate scaling of the amplitudes. We therefore consider now the solutions of 

dfUJ + divj. (wu) = 0, 
1 ^ 

h':it)^[h's)^ik{t,h.,it))y, 

3.4 A regularized version of the system 

Such an issue is quite similar to the one which consists to obtain the Vlasov-Poisson system as a mean-field limit 
of Newton's equations for charged particles. This last problem is still open at the time of writing0 but some 
results are available in the simpler setting where H is assumed to be W^'°°, in the spirit of the famous paper [TS] 
by Dobrushin (let us also mention here Braun and Hepp [8] and Neunzert [27]). We will therefore consider the 
equations (P)-© where u is given by 

u := K[cj + p] and p -.^ [ fd^ (30) 

where K is defined by 

K[g]{x):= [ H{x-y)g{y)dy, 
with H is in W^^°°{R'^) and satisfies H{0) = 0. 

Theorem 1. Assume H g W^'^^iM?), satisfying H{0) = and K defined as above. We have the following. 

(a) Assume that (wo, /o) A^i(M^) x 7^i(R^ x M^). Then there exists only one weakly continuous curve {uJt, ft) 
in "^"([0,00); A<i(R2) x Vii'^^ x M?)) solution of the system ^-^-^ with {ujq, fo) as initial data. 

(h) Moreover we have the following stability property. Consider two solutions pi := (cji,/i) and p2 {^2, f2) of 
the system ^-^-^^ associated to two initial data := (wqi/o) and pi := (wg,/o) m7Wi(M^)xPi(M^xR^). 
Then, for any t ^ 0, 

WMt),p2{t))^e^'''W,{pl,pl), (31) 

where C > depends only on ||-ff||Lip and |aJo|(K^). 

(c) Finally if we assume that (aJo,/o) *s also in Lip(R^) x Lip(M^ x R^); then the corresponding solution {ujt,ft) is 
m Li^^([0,oo);Lip(R2) x Lip(R2 x R^)) 

Above Wi is the Wasserstein distance defined in Dcfinition[5] In the proof of Theorem[T]wc will use the following 
notation: when ly is in A^i(R^), we denote Ulv] the vector field defined for (a;,0 G R^ x R^ by 

UM{x,0 (e,(e-A'M(.^))^) eR' xR2. 

We will use the two following lemmata. 

^In the case of the Vlasov-Poisson system, Hauray and Jabin have recently succeed to improve the previously known results about 
the mean-field limit to some cases where the interaction kernel can be singular cf. |20) . Yet their approach docs not cover the case of 
the Newtonian kernel appearing in JSj. 
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Lemma 4. Let v he in A^i(R^). Then the vector fields K[iy] and U[iy] are Lipschitz and 
\\K[iy]\\up ^ Wm^)\\H\\up and||C/H||Lip^ |i^|(M2)max(||i/||Lip,l). 
Moreover Kliy] is uniformly bounded 

Wm')\\H\\L^- (32) 
Proof. This follows easily from convolution properties, referring to the definition of K[-] above. □ 
Lemma 5. Let vi and 1^2 be in Mi{M.'^). Then 

lit/N-C/HIU- = ||A'M-if[t^2]||L- < \\H\\upW,{iyui^2)- 

Proof. We have 

uWi] - c/H = (0, -{KWi - i^2]{x))^), 



so that 



|C/H-C/h]| = \K[iyi~i^2]\ 

= \H^{vi-V2) 



\H\\up 

It then remains to use Proposition [T] to conclude the proof of the Lemma. □ 

Proof of TheoremU[ Denote no (wc/o) G MiiM."^) x ■Pi{R'^ x M^)^ 

Consider the subset C A^i(R^) x x M^) of elements (w, /) such as a; is compatible with ujq. Then, for 

T>0,T:= 'S'°{[0, T];To) endowed with 

Wl{^ll,^l2)■■^ sup Wi{ni{t),fi2{t)) (33) 

te[o,T] 

is a complete metric space. We first prove the local existence and uniqueness of solutions to the equations ([T])-©- 
pop by using the Picard-Banach theorem on J- with T > small enough for the mapping 7^^, we are now going to 
describe to be a contraction for the previous distance. 

The mapping 7^^ is defined as follows. Denote by tt^ : x — >■ R^ the projection on the first variable. 
Given /xq := (wo,/o) in -^0 and /i := (w,/) in T, we define p[f](t) := 7r^/(i) = (fit), 1)^, so that uj + p[f] is in 
'r°([0, T]; 7Wi(R2)). Then thanks to LemmaH the vector fields K[uj + p[f]] and U[uj + p[f]] are Lipschitz, uniformly 
in time; and we can therefore define 

• a unique continuous flow map 0^ deflned for {t, x) G [0, T] x R^, with values in R^, given by 

0^'(x) X + f K[uj{s) + p[f]{s)\ {^^{x))ds, (34) 
Jo 

• a unique continuous flow map defined for (i,x,^) G [0,r] x R^ x R^, with values in R^ x R-^, given by 

E'^{x,0:={x,0+ f U[uj{s) + p[f]{s)]{j:^{x,0)ds. (35) 
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Moreover, for all t S [0,T], /i^ and are Lipschitz and we have the following estimate 

ll^niLip + llSfllLip ^ cxpj^ max(||H||Lip,l) [|c^(s)|(M2) + l]ds[> (36) 
= exp 1^ max(||i?||Lip, 1) [|wo|(K') + l]ds 
We then define Tf^a[^] as a function defined on [0,T], with value in 7M(R^) x by 

Actually, 7^^, takes its values in J". 

Indeed, since pushforward transformations conserve the mass and transport the Jordan-Hahn decomposition, for 



anyie [0,T], ^ (t)) e ( 



'), u'^it) being compatible with wq. Moreover using p2p we have that 



for any t G [0,r], \\K[u} + p[f]]it, ■)\\l'-'{r^) < " with a ||iJ||ioo(R2) (1 + |wo|(M^)), so that we see from ([Mll-dSa) 
that \(j)'^{t,x)\ ^ |x| + at and |S^(t,x,^)| < (|a;| + |^|) + t{a + 2|^|). Since ujq and /o arc finite measures having 
finite first moments, the previous inequalities insure that so do and f^{t), for all t e [0,r]. Finally in order 

to prove that T^,, [^] is in J-' it only remains to stress that the continuity in time of i/)^ and yields the continuity 
in time of T^j, [/x] . 

Let us now observe that a fixed point of the mapping corresponds to a weak solution of the equations 
([!])- with IJ.Q = (ojo, /o) as initial data. Now let us see that 7^,-, is a contraction for T small enough; consider 
yUi := and /i2 := (w2,/2) m J". 

Let us define, for t in [0,r], 



m 



10; 



|c?|wo| and X{t) := 



We have, from 



X\t) ^ 



thanks to 



Now we use the triangle inequality to get X'{t) ^ Xa{t) + Xb{t) with 

A,(i) := / K[uj2{t) + p[f2]{t)] o ^t' - K[u2{t) + p[/2]W] ° d\u:^\ 



We have, for any t in [0, T], 

Wi[u:i{t) + p[h][t),u:2{t) + p[h][t)) ^ Wi{ooi[t),u:2[t)) + Wi[p[h]{t),p[f2]{i)) 

^ Wi{pi{t),pL2{t)), 

where we used Lemma [T] for the first inequality and Lemma |3] for the second one, since p[fi\{t) = T^^fi{t), for 
1 = 1,2. Thanks to Lemma [5] we therefore obtain that 



Xait) 



K[uj,{t) + p[h]{t)\ - K[u2{t) + p[/2](t)] 



|i^||LipW^l(MlW,M2W)ko|(K')- 



On the other hand using Lemma HI we obtain 

Xb{t) < ||iJ||Lip(l + |^o|(M2))A(f). 
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From the equations ([34 |) - ([35l) we infer that A(0) = 0. Then the Gronwall lemma leads to 

A(t) ^ C / e^(*-^'l^i(Mi(s),/i2(s))ds, 



for any t in [0,T], where C > depends only on ||-ff||Lip and |cl>o|(IR^)- Since, for any couple {i'i,i'2) G A^i(M^) x 
A^i(M^), |C/[i^i] — [^[1^2] I = the same computation gives a similar bound for A(t). 

Now, for any t in [0,T], denoting by the index ip a supremum over the unit ball of Lip(M^) and the index Tp a 
supremum over the unit ball of Lip(]R^ x R-^). 



= sup / (^d((/)fla;o - (^flwo) + sup / d(Sj'it/o - Sfi/o), 

= sup / (</jo(/)^i -(^o0f )dwo + sup / (V'oEj'^ -V'oSDd/o, 

A(i) + A(t) 

^ C [ e^(*-^)M^i(Aii(5),/.2(s))d5, 
Jo 

Therefore 

Wi(r[Aii],rM) ^ (e^^-i)Wi(Mi,M2), 

so that 7~ is a contraction for T small enough. This smallness condition only depends on the total variation of lhq 
and on the Lipschitz modulus of H. Since uj{t) remains compatible with wq at all time, it has in particular always 
the same total variation. We hence infer the global existence part of the result by an iteration process. 

Let us now prove the stability estimate (PT|) . So let us consider two solutions fj,i := (wi, /i) and ^2 := (^2, /2) 
of the systems pl-(P|)- ([5n|) associated to two initial data /ij := (wo,/o) and /Zq := {ujq^/q). Then, we have the 
following where t in [0, T] is understood, 

wMt),fi,{t)) = w^i(</>rH,</>r^^o')+w^i(sr^/o,sr^/o)- (37) 

Let us first focus our attention on the first term in the right hand side. We have by the triangle inequality that 

Now the first term can be tackled as previously whereas the second one can be bounded by using Lemma [5] and 
the estimate We proceed in the some way for the second term in the right hand side of ((57)) . This proves the 
desired estimate. 

The last statement of the Theorem follows classically from basic transport theory. □ 
3.5 Mean-field limit for regularized kernels 

We infer from Theorem [1] the following result about the mean- field limit in the regularized case where (|30p is taken 
instead of ([3]). 

Corollary 1. Assume that (wo,/o) is in Mi{R'^) x Vi(M? x R"^). Let {h°,h}),en* e (R^ x R^)^' such that the 
sequence of empirical measure (/o^)jV6N*7 defined by 



1 " 



i=l 

satisfies Wi{f^ , fo) — > when N — > +oo. Let us denote by {'^^ , f^)N£W* , ("^i /) solutions 

fN 

lo 



respectively associated to the initial data (wq, /(j^)ArgN* , (i^Oi/o) given by Theorem\^ Then 
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1. for any t > Q, for any N ^ \, 



where 



N 

4=1 



N 

\hi,N{t),h'. „(*)), 



KAt) = [K^j,{t)~u''{tM.N{t))Y. (38) 

1 ^ ^ 

K[u^''] + j.Y.H^--^Mt))^ (39) 



;0 . 1\ 



N 

N 

{h,MO),KNm = {hlh\). (40) 
2. for any T > 0, yVi(/iAr,/i) — > u;/ien — > +oo, where Wi is the distance defined in p3p . 

Let us recall here that (0) = so that the velocity um can be seen as a regularization of the velocity u of 
Section (|3.3p and of the velocities Vi as well. Hence the equations ((38)) - (|40| . together with the following transport 
equation for the vorticity: 

9(0; AT + div(a;Ar ujy) = 0, 
can be seen as a regularization of the equations of Section p.Sp . 

It is also possible to obtain the following result for the case where the fluid vorticity is also discretized. 

Corollary 2. Assume that {uJa, fo) is in Mi{M?) x "PiCM^ x M?). Let {x° , a,, h° , hl),(zw* € (M^ x M x x m2^n* 
and consider the sequence (w^, /(j^)jv6N' € A^i(R^) x ViiM."^ x K^) defined by 



UJ 



1 " 

i=l 

where {i G : ±ai > 0} and jS^ := J^iei^'^i- Assume that Wi{u!q ,ujo) + Wi(/(j^,/o) — > when 

N — > +00. Lei us denote by '■= {uj^ , f^)N£fi' , l-J- '■= {^,f) the solutions respectively associated to the initial 
data {lJq , fi^ ) N ^fi* , (wqj/o) given by Theorem[J\ Then 

1. for any t > Q, for any N ^ 1, 



W = r,+ «^'5^..iv(t) + o- 2^ a«o^.,«(t) <'^d fN{t) = — ^5(h^ ,,(t),h[ „(t)) 



where 



4,Nit) = u''it,X^,Nit)), (41) 

KN(t) = (/<w(t)-Sf^(t,/i,,w(i)))^, (42) 

^I^ = K[u;^] + ^ E Hi-~h,Mt)), (43) 

(a;,,jv(0),/i.,w(0),/iU(0)) = {xlhlhl). (44) 

^. /or any T > 0, yVi(/iAr,//) — > when N — > +oo, where Wi is t/ie distance defined in p3p . 

The system (|4T |) - (|44| describe the dynamics of 2A^ pointwise vortices localized in {xi^N{t))i!^i^N and in 
{hi^N{t))i^i^N under pairwise regular interaction, the first N ones being massless. 
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3.6 Stability of the hydrodynamic regime 

Given a map v : M'* — > R**, one can define its "graph" function 

t" ■.W' — > X 

X I > (.T, v{x)). 

For any measure p e A^(R''), the pushforward measure r^p is hence a weU-defined measure on R*^ x R''. In 
the particular case where u is a continuous function, one may write, for any bounded continuous test function 



(r|p, (/)) = (p, o r") = / (j){x,v{x))dpix), 

SO that the measure r^p can hence be seen heuristically as the product p{x)Sy(^,j.-j (^) and may be called monokinetic. 
The generalization to the case of time-dependent measures/functions is straightforward. Actually, another corollary 
of Theorem [1] deals with the case of solutions taking the previous form. Formally, in a time-dependent frame, if 
/ := T^p the equations H])-© reduce to the following system for the unknowns uj{t,x), p{t,x) and v{t,x): 

dfU! + diYx{uju) ~ 0, 

dtp + divx{pv) = 0, 
dt{pv) + divxipv 'Si v) = p{v~u)^, 

where u is still given by 

u = K[uj + p]. (45) 
When p does not vanish the third equation in the system above can be simplified into 

Here again we will deal with the case of a regularized kernel substituting the law 

u = K[uj + p] (46) 
to (j45p . For such a system we have the following result of loeal-in-timc result. 

Proposition 2. Let {ujo,pq,vo) G 7Wi(R^) x ViiM."^) x Lip(R^). Then there exists T > and {u,p,v) e 
L°°(0,r;7Wi(R2) X Vi{R^) x Lip(R2)) solution of 

dtto -\- A\v x(y>u) = 0, 
dtp + diYxipv) = 0, 
dtV + v-VxV = {v~u)^, 

with (w(0), /9(0), w(0)) = (wojPoi'^^o) and u given by (|46p . 

Proof. The proof is an easy application of the method of characteristics. □ 
Then we have the following. 

Corollary 3. Assume that (wo,Po,wo) e Mi{M.'^)x'Pi{M.'^)xUp{M.'^) . LetT>0 and {oj,p,v) e L°°{0,T;Mi{M.'^)x 
7^1 (M^) X Lip(R2)) be given by Theorem\M Consider fo := r^o^po G Pi(R^ x R^). Let be given a sequence (/,f )Ar^i 
is in ViiR"^ x M^j s^^/i that Wi{f^ Jo) when N +oo. Let us denote by p.^ := {uj^ J^)n^i, p := (w,/) 
the solutions respectively associated to the initial data {ujq, f^)i\[^i, (a;o,/o) given by Theorem\^ Then, f actually 
equals to T^p and, for the time T given above and the distance Wi defined by (j33p . yVi{p^ , /i) — > when N — > +oo. 
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4 Existence results 

4.1 Weak solutions 

The main resuh of this section is the fohowing 

Theorem 2. If loq is in is in {L°° C\L^){^ x M?) such that the kinetic energy of the dispersed 

phase is finite: 

I fo{x,m?dxd^<+^, 
then for any T > there exists at least one weak solution 

pS[l,oo[ 

to the equations ((T])-([3]). Morover the kinetic energy of the dispersed phase of this solution is finite at any time: 



vte[o,r], / f{t,x,0\i\'- 

JR2xR2 



dxd£, < +00. 



Proof. Let us introduce some notations. Given two functions g(^) and and a nonnegative number a, we 

denote: 



maig):= giOrn^d^ ixnd AUh) := hix^Ol^^dxd^ 
We have the fohowing useful result 

Lemma 6. If g G L°°(R+ x x M^) (positive) such as m~f(g)(t,x) is finite almost everywhere in time-space, then 
for all a < J we have, almost everywhere in time-space 

ma,{g){t,x) ^ C{\\g\\L^^ + l)m^{g){t,x)^ . 

Proof For all i? > 

m^{9){t,x)= [ g{t,x,mrdi+ I g{t,x,OWd^ 

^C\\gh^R'^+' + j^m^ig)it,x), 

and the real R = m^(t,x)^i+^ does the trick. □ 
Now consider (ct'Q,/g'') a smooth compactly supported approximation of (a;o,/o)- Recall the definition 

^ ^ 27r |a;|2 

Consider a regularizing kernel (<Pn)neN of pair functions, in such a way that for all n € N 

Hn := H^ifn e VK^-°°(M^), and H„(0) = 0. 

Theorem [1] allows us to construct for every n G N a global strong solution (w", /") to the system 

dtuj'' + u" • V^w" = 0, y{t, x) e [0, T] X (47) 
uj''{0,x) =u^{x), VxgK^, (48) 

9tr + e-v,r + (e-ii")^-ver = 0, y{t,x,o&[o,T]xR^ xr^, (49) 

f^{Q,x,0^f!>'{x,0, V(a;,e)eR2xM2^ (50) 
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where m" := Knlp^ + w"] and Kn is defined by Kn[g\ = K[g\ -kipn = Hn with ((^„)„gN- Note that w" and /" 
are compactly supported (and the later positive), and we have, for all n G N, 



||r||L~ = ||/o"l|L- S^ll/olk-, 



(51) 
(52) 

(53) 



From now on. will denote a positive constant (that may vary from line to line) depending only on the initial 
data and T. 

Denoting :~ / /"d^, j" := / /"C*^^, and applying lemma [5] with respectively a = = 2 — 7 and a = 1 = 7 — 1 
we get 



||p"(t)||i. = l|mo(r)(t)|U2 CinA^2(^)(^)^/^ 

i1j"WIIl4/3 iimi(r)(i)iiLV3 < a„M2(.r)(t)3/4. 

Using Holder inequality and ([551) - ([5^ we get 

Thanks to the Hardy-Littlewood-Sobolev inequality we get 

||«"(t)|U4 <;qip"(t)+.."(t)|l^.,3 «;QTJ1 + M2(r)(i)'/'] 

Now, multiplying the kinetic equation by and integrating in space- velocity gives 



(54) 
(55) 



(56) 



dt 



M2{rm 



r{t,x) ■ u"^{t,x)dx < 2||j"(t)|L4/3h"(i)|| 



and hence, using (|55p . 



so we get 



4 

dt 



\u-{s)\\Us 



and using ((56)) we have finally 



Ahirm^ci 1+ / M2{n{s)ds 



which, via a linear Gronwall lemma, allows us to conclude that (Af2(/"))„gN is bounded in L°°{[0,T]). Using 
((55)) and ((5^ we hence deduce the boundedness of {p")neN, (i'OneN and {u^)neN respectively in L^(L^), 
LriL'J') and Lr(i^)- 

Now (|47p gives dtuj" = — div2;(u"a;"), and because (a;")„gN is bounded in L'^{Lt^^), we have that (9ta;")„gN 
is bounded in L'^{'T^^{M.'^y) ^ ^t°{H~\^^) for m large enough. Analogously, we have 9tp" = —divxj" and hence 
(atp")„6N is bounded in the same space L^iH'l'^J. Now • * V^G sends iJ~™(M2) in another iJi;^«(M2)^ so 
that we have finally (9tw"),igN is bounded in L'^{H~l^^. But, from Calderon-Zygmund theory we know that 
11 VM"(t)|L4/3 ^ C\\p^{t) + a;"(t)j|^4/3 ^ Cj5^, and since (w")„gN is bounded in Lf'{L^), it is also bounded in 
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L^{W^\J^ ). Wc hence apply Aubin's lemma to get (up to an extraction) the strong convergence of (u"')„gn in 
L^(L^ jqj,) for all p < 4. At this stage, we cannot use a weak/strong type argument to pass to the limit in both 
equation: the bounds (j51[) and (|52[) insures that up to an extraction (/")„gN and are weakly converging 

respectively to some / and uj in L'^{L^^) — ★ and L^(Lf/'^) — but (M")„gN does not a priori converge in L]{L'^) 
so that we cannot directly pass to the limit in the nonlinear term of (I47p . However it follows from the stability 
part of [12] that (a;")„gN and (/")„eN respectively converge in ^t^^l^Lt^^) and ^t{L^ for aU p G [1, oo[. We hence 
pass to the limit in both equations and obtain by the same time a weak formulation including this time 

the initial data. This concludes the proof. □ 



4.2 Classical solutions 

In the spirit of }14| we construct regular solutions to our system, taking some regular initial data, with a kinetic 
phase having sufficient decreasing at infinity in the velocity variable. More precisely we have the following result: 

Theorem 3. Assume ujq € ^^^'^(M^) n W^'^^iM."^) and /o G ^^^'^(M^ x R^) j^thg latter being positive), such as 

(1 + ler^HlM + \Dfo\) e i°°(R' X m"), (57) 

for some real number 7 > 2. Let T > 0. Then there exists a classical solution 

{ljJ) <E i°°([0,r];M/i^i(M2)) X L°° {[0,T];W^'\K.^ x R^)) 

to the equations satisfying furthermore 

(1 + \er^H\fo\ + \Dfo\) e L-([0,r];L-(R2 x R^)), 

Remark 3. Taking W^'™" initial data and replacing ()57p by a similar assumption taking into account all the m-th 
first derivatives lead to even more regular solution. This is done in [14^ , the proof being quite close to the case 
m = I. 

As wc follow closely the method of [14] we will only provide a sketch of proof. More precisely we will only prove 
some appropriate a priori bounds for the approximations given, as in the previous section, by ()47p - (|50p . 

Proof. We consider again the solutions of (|T7)) - ([5n]) and we define as in [T3] 
we have, using (|49)) 

atF" + ^ • v,r" + (c - u")^ ■ v^r" = 7^"^ • e(i + lei')'''"'^^'/", 

which is just a classical transport (Vlasov) equation with a source term. We hence have, thanks to the assumption 
on the initial data 

\\Y-it)u^ < 11/0(1 + \er^'\\L^+i f 11(1 + \ep-''>'^ns)\\LAu''{s)\\L^ds. 

Jo 

But It" := Kn[p" + cj"] so that 

||u"(s)|U. ^q|p"(,s) + c.(,s)||^/^!|p"(,s) + c."(s)||ie 

^ [IIpoIIli + \\u:o\\L^f'[\\p"is)h^ + \\uo\\l^Y^\ 

where we used for the first inequality a classical interpolation estimate true for the Biot-Savart operator and hence 
for its approximation /\„, uniformly in n. The second inequality is a direct consequence of equations (|47p and (|49p 
and the definition of fg and ujq as rcgularization of the initial data. We hence use, as in [14] (see appendix) the 
interpolation estimate 

IIp"(.)IIl~ c,\\f-is)\\^;!-'^/''\\Y-is)f^2. 
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We have also 

11(1 + i^n(''-i)/V"(s)iu- ^ c,\\ris)\\]^2\\Y^\s)\f£-'^/\ 

so that we have finally a Gronwall estimate 



1 



f\\Y-{s)U^ds 
Jq 



We have hence at this stage bounds for (F")„, (/9")„ and (u")„. Now introduce, again as in 

z"{t,x,o (i + i^i'r/'v,,^r(i,.T,o. 

We have the vectorial equality 
where A" G is the matrix 

O3 V,u"-^ 



-Ids -V^e' 
and hence, 

Using first a logarithmic estimate a la Beale-Kato-Majda, and then the previous bounds, we get: 

||V,u"||l=o < c[i + + c."Gs)||i»(l + In {1 + ||V,p"(s) + V,L."(s)|U=e }) + + c."(s)lh 

<C[l+ln{l + ||V,p"(s)||Loo + ||V,c."(s)|U=o}] 

In the same way 
verifies 

so that we get eventually for some constant C, letting Q;„(t) ;= 



Oln{t)^C 1+ / Q!n(s){l + lnQ!„(s)}(is 



(58) 



VK"(t)|| 



which classically implies a bound for (a„)„ in L^, and hence a bound in L'^{L'^^ for the sequence of matrices 
(A")„. Integrating ([58]) and using the last bound we have also 



iv,,^r(s)iiiids 



that is (V:e,j/")„ bounded in Lf{Ll^^). From all the previous bounds we get (up to several extractions) convergence 
in — w-k frame. From this the proof of Theorem [3] follows exactly the same lines than Theorem 1.1 in [T!]. □ 
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5 A uniqueness result 

Let us give here a uniqueness result which extends the result |31| by Yudovich about the incompressible Euler 
system and the result [23] by Loeper about the Vlasov-Poisson system. 

Theorem 4. If ujq is in {L°° n L-^)(R^) and Jq is in A^J(M^ x M^) then, given T > there exists at most one 
weak solution 

{uj, f) e '^^°([0,T]; n l1)(M2) _ y; ^ xM + iR^ x R^)^ 
to the equations (P-© such that p is in L°°{{0,T) x M^). 

Note that Theorem S] applies to the case of mono-kinetic solutions, that is for a density / of the form given 
in Section 13.61 The proof of Theorem U] follows closely the method of [33] . Yet we provide a proof for sake of 
completeness. 

Proof. Let be given two solutions and (w2,/2) of the equations (H])-® such that pi := mo(/i) and p2 := 

"^o(/2) are in L°°{{0, T) x R^). Accordingly we define Ui := K[uji + pi], for i = 1, 2. 

Recall the following classical result (sec for example [H])- 

Proposition 3. Let g £ L°°(M^)ni^(M^). Then K[g] is bounded and Log-Lipschitz on M^, that is there exists C > 
such that \\K[g]\\L^ ^ C, and such that for any {x,y) € x with \x — y\ ^ i, there holds \K[g]{x) — K[g]{y)\ ^ 
-C\x - y\\n\x - y\. 

It therefore follows from the assumptions of Theorem |4] that the vector fields ui and U2 are bounded and 
Log-Lipschitz on M^, uniformly in time, and so are the two vector fields, for i = 1, 2 

: [0, T] X X m2 — ^ jg2 ^ ^2 

^ (^^,{(^u,{t,x))^y 

Because of the mentioned regularity, referring again to jl2| . we can define, for i — 1,2, two unique continuous flow 
map (f)' : [0, T] x and : [0, T] x x ]R2 ^ r2 x such that 

4)].{x) := <t)\t,x) = X + [ Ui{s,(j)'{s,x))ds (59) 
Jq 

mx, ■■= ^\t, X, = {x, 0+1 U\s, ^\s, X, 0)ds. (60) 

JO 

For i = 1,2, define two vector fields X\^\ defined on [0,T] x x ^j^-j^ values in such as E'' = (X*,S*) 
that is to say that X'^ and are respectively the characteristic curves in the space and velocity variable. Denote 
as before Xl{x,Cj := X''{t,x,£,) and Sj(a;,C) := 'E}{t,x,Cj for all {t,x,Cj G [0,T] x x M^. 

We use the following result cf. Theorem 1.2 in [5] or [T] for a rather different proof. 

Theorem 5. Let V{t, x) : [0, T] xM'' he a Borel vector field such that there exists C S L^{Q, T) such that for all 

t in [0,r], !|V(t)||Loo <; C{t), and for any {x,y) eR^xW^ with \x-y\ ^ \, \V{t,x)-V{t,y)\ < -C{t)\x~y\\n\x-y\. 
Let i^iQ G jM(M'^). Then there exists only one weakly continuous curve fit in '^°([0, T]; A^(R'') — w*) solution of 
dtpt + div2;(/UtV) = 0. This solution is given by fit = ^Y#f^o, where is given by flow associated to V at time t 

by the formula <^>^ {x) := $^(t, x) = .t + ( V{s,<^'^ {s,x))ds. 

Jo 

As a consequence, there holds, for i ~ 1,2, t e [0, T], 



{^i, fi, Pi){t) = (0J#wo,Sj_^/o,^t#/o), (61) 

i li 



u;±,k|)(i) = (0j#wo±,0j„|wo|), (62) 
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so that at any time t £ [0,T], 0Ji{t), fi{t) and pi{t) are respectively compatible with uj2{t), f2{t) and P2{t). For 
two time-dependent measure (always compatible) we will often use the shorthand 1^2(1^1, '^2){t) ■= W2{i'i{t), i'2{t)). 
Finally we introduce the quantities 

Qit) -.^ I |Ei - S^pd/o, Q{t) / - ct>l?d\uj^\. 
Let us now recall the following result cf. [53]: 

Lemma 7. Let ijji : M'' — > R'', for i ~ 1,2, be two homeomorphisms . Let v E A4 + {M.'^) and denote, for i ~ 1,2, 
i^i '■= Then, for vi and V2 are compatible and 

Thanks to this lemma we obtain, for any t £ [0, T], 

VF2(/i,/2)(i)' Q{t) and W2iui,uj2){tf + I^aK", ^2") W ^ QW, (63) 

so that it suffices to prove that Q{t) and Q{t) vanish to get that (ui, fi) = [102, f2)- 
Differentiating in time we get 

Q'{t) - 2 / (Si - E?) • (9il]i - dtT?,) dfo, 
Q'{t) - 2 / ~ 4,1) ■ {dtc^l ~ M't) d\u;o\. 

JR2 

Using now ([60]) and ([59]) 

Q'it) = 2 [ {Xl~X^)-{^\~^1)dfo 

-2 f (El - S2) . {u,{t,Xl) - U2{t,X^))^dfo, 
Q'{t) = 2 f (01 - </)?) . {ui{t, </)!) - U2{t, 0?))d|wo|. 

JR2 

Using the Cauchy-Schwarz inequality, we have that the three terms in the right hand sides above are respectively 
bounded by Q{t), by 

2Qit)U [ \m{t,Xl)~U2it,X^)\^dfoy, (64) 

and by 

r - 
2Q{t)U \uiit,4)l)-U2it,4)^t)\'d\iOo\y. (65) 

^ JR2 ^ 

Let us introduce the quantities 



Tiit) := 


/ \u2{t,Xl)~U2{t,X^)\^dfo 

jR2xR2 


T2{t) := 


^Ju2{t,X^)-U,{t,X^)\^dfo 


Ti{t) := 


[ \u2{t,cl)l)-U2it,4^^t)fd\uJo\, 

JR2 


T2{t) := 


/ \u2{t,C^l)~Ui{t,(f>l)\^d\uJo\. 

JR2 
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Then the terms in (|64l) and (|65|) are respectively bounded by 

2Q{t)HTi{t)^ +T2{t)^) and 2Q(t)^ + TzW^)- 
Using (|5T|) and (p^ we have 

T2(t) = / |u2(t)-Hi(t)|'dpi(t), 

Mt) = [ \u2it)-Uiit)\^d\cJi{t)\. 

The crucial estimate in Loeper's method is the following. 

Lemma 8 ([23], Th. 4.5). Consider two compatible signed measures on R^, and ^2, both having bounded 
densities with respect to he Lebesgue measure on M^. Then 

\\K[fii] - K[fi2]\\L'^ ^ (max{||^j+||L^, ||/i+||Loo})5 W2ifit , t^t) + (max{||^j"||Loo, W2(p^,ti2)- 

In particular we get here, for every t G [0, T] 

- U2mL- ^ C[W2{iO+,U+) + W2{lO^,L02) + ^^2(^1,^2)] W, 

where the constant C depends on max(||cji||ioo, ||Lj2i|L°°, ||/52|1l=°)- Modifying C into another constant 

depending only on the same quantity, we get: 

max(T2(i),T2(0) ^ c[w2{uj+ + W2{u^ + W2{pi, P2)\\t). 

Using now (j63p we obtain that for an appropriated constant C, 

max(T2(i),r2(t)) ^ C{Q{t)+Q{t)). 

Let us now tackle the terms Ti and Ti. From the L°° bounds on ui, U2, and the equations (j59p and (pO]) we infer 
that there exists > (depending only on the initial data through their norms in n L°°) such that 

110"^ - 0^llL°°([O,t*]xE2) + - -'^^||L~([o,t,]xE2xE2) < -. 

We may even assume that the same estimate holds for Q + Q on the same interval. Then, since the vector fields 
Ui and U2 are Log-Lipschitz on K^, uniformly in time, we have, for all t e [0,i^,], noting g : .t i— > a;(hia;)^ 



Tiit) ^ c[ \ixl~xl)Hxl~xf)\'dfo^^ [ g{\xl-xf\ 

jR2xR2 4 jR2xR2 



jR2xR2 

c 



T,{t) ^ C m-4>^,)H<l>\~4>l)?d\u^\^- g{\4:\-4>l?)d\oo^\ 

JR2 4 Jr2 

Since g is concave for ^ x ^ i it follows from Jensen's inequality that (after renormalization) 



2|2 f^/o 



l/olL 



Since g is increasing for ^ x ^ i we get on [0, t^\ 
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Thus, doing the same for Q, we get eventually 

Ti ^ C/„0(l-lnO)2 
Ti s=: a„Q(l-lng)2 

Collecting all the previous bounds we obtain an estimate of the type, on [0, t*] 

g' + O' ^ Cf,,^, [Q + g-QlnQ-Qlng], 

and by convexity we get finally 

Q' + Q' ^ Cf„,^, [Q + Q-iQ + Q) HQ + Q)] , 

Since Q{0) = Q{0) = 0, this allows to conclude that Q and Q vanish on [0,t^] and thus on [0,r] by iterating the 
above argument. □ 

6 Poisson structure 

The goal of this section is to exhibit the Hamiltonian structure of the equations (d))-©, following [3]. We will be 
quite formal here, leaving aside the regularity issues. One therefore will think at some functions (w, /) which are 
smooth in t, x, ^ with a nice decreasing at infinity. First we endow the manifold V of the pairs (a;, /) with a Poisson 
structure, that is to say, we endow P with a bracket {•,•} acting on 't^°° functionals F : V ^ M., bilinear and 
skew-symmetric, satisfying the Jacobi and the Leibniz identities. Here this is obtained by setting, for any smooth 
functionals F,G on P, 

{F, G}:^ f f W^Ff ■ WfGf + f f {Ff, G/},,^ + f w ^ .,F.^ ■ W^G^, 

JR2xR2 jR2xR2 JR2 

where F^^j and Ff denote the gradients with respect to lu and / of a functional F in , given by the formula 

f)h . hm Pi- + ehJ)-Fi.J) f ^^^^^ ^~ ^ f + sH) - /) 



and the bracket {•, ■}x,(, stands for 

V,/- V^g- V^/- V,g. 

The properties cited above are clear from this definition. 
Let us now define % by 

2H^ I efit,x,Odxd^- [ Gix-y){uj + p){t,x){u; + p)it,y)dxdy, with p := / M, 

jR2xR2 jR2xR2 JR2 

which can also be seen as a functional on the manifold V, and this endows the system with a Hamiltonian structure 
in the following sense. 

Proposition 4. When {uj,f) solves the equations ([I])-© we have for any smooth functional F, the ordinary 
differential equation 

^F = {F,n}, (66) 
where F and H in stand respectively for F(ui, f) and H^uj, f). 
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Proof. According to the chain rule, we have 



dt ^ ./b2 dt ./b2vB2 dt ^ 



Using the equations (HJ-© we arrive to 
d 



On the other side the derivatives of H are 

n^^-[ G{--y){u: + p){y)dy, = ^ - [ G{- ~ y) {lo + p){y) dy, 

JR2 ^ JR2 

Then, integrating by parts, we find 

JR2 JR2 

since u = V-^Huj- 

On the other side, since ^ = V^H / we have 

h= I fV,FfV^Hfandh = f fV^FfV^Hf- f fV^HfV^Ff, 

JR2 X R2 JR2 X R2 JR2 X R2 

what yields ([66|). □ 

As a simple corollary of (|66p and of the skew-symmetry of the bracket we get that 7i is conserved by the 
solutions of dU-®. 

7 Massless limit 

In this section we investigate the behavior, when e — > 0+ of the system 

dtuj^ + divx{uj'u^) = 0, (67) 

dtr+dw,{ro + -divdn^~^')^) = o, (68) 

e 

where 

:= K[uj' + p'] and p" / f'dS,. (69) 

JR2 

This corresponds to the equations (II])-®-® with an extra factor which encodes a regime where the particles 
immersed into the fluid are lighter and lighter as e — > 0"^. 

We are going to prove that, in this limit, the equations (|67p - (p5|) - (|5ni) degenerates into the incompressible Euler 
equation: 

dtu + dwx{u®u) + \/p = 0, (70) 
div^u = 0. (71) 

An heuristic way to guess the previous limit is the following: we infer from the equation (|68p that in the limit 
£ — > 0+ the density of particles becomes monokinetic with a velocity ^ = u so that 

f ■■= I r^di (72) 
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converges to pu, where p and u stand for the respective hmits of p^ and . Therefore the right hand side of the 
equation ^ vanishes, again formally, what yields ((7n)) - ([7T|) . 

This problem is very close to the so-called gyrokinctic limit considered in ^ with here an additional coupling 
to an Eulcr-typc equation. We will actually follow the strategy of based on the use of a modulated energy, with 
a few modifications. In particular we will make use of the notion of dissipative solutions of the Euler equations. 
This notion was introduced by P.-L. Lions in [52] and used in [5] with a slight modification in the definition. Here 
we will perform another slight modification by extending Brenier's definition to the case of infinite energy, since 
arguably in 2d the case of a incompressible perfect fluid of finite energy is too restrictive. For instance, it is easy 
to see from the definitions dU and ^ that, given a smooth compactly supported function g from to R, we have 
that 



K[g\ is in if and only if / g{x)da 

JR2 



0. (73) 



To deal with vorticities which do not satisfy this last condition, we follow [H] (see also [23]) by defining the following. 

p+oo 

Let a € R. There exists ga G i^(]0,oo[;R) such that 2tt / gair)rdr = a. We then define Ha ■= K[ga{\\ ■ ||)] 

Jo 

and the space 

Ea Ha + L^{R ), 

where L^(R^) denotes the divergence free vector fields in i^(R^). Observe that Ha is a smooth stationary solution 
of the the 2d incompressible Euler equation (| 70 1) - ([71]) . Finally the finite energy case corresponds to a = : thanks 
to ([72]), there holds Eq = L'^. 

Definition 6. Let a S R and uq € Ea- We say that a divergence free vector field u e L°°((0,T); Ea) is a 
dissipative solution of the incompressible Euler equation (|70p - (l71l) with uq as initial data if for any smooth vector 
field V e '^°{[0,T];Ea), such as A{v) ;= dtv + v-Vv e L^{{0,T); L^{M.^)), and d{v) £ L^{{0,Ty, L^{«.^)) (where 
d{v) is the symmetric par^ of D{v)), for almost every t G [0,T], 

\u{t,x)~v{t,x)\^dx s^, [ \uoix)-v{0,x)\^dxe^p [ 2\\div{0))\\de (74) 







2 / / A{v){s,x){v~u){s,x)e^p{ / 2\\d{v {d))\\d9 } dxds, 



where \\d{v{9))\\ is the supremum in x of the spectral radius of d{v){9,x). 

Remark 4. Just as in ]2'^ . let us note that by a standard regularization procedure it is sufficient to verify estimate 
(|74p for (only) all smooth (in time/space) vector fields v satisfying v{t) e Ea and curlu(t) is compactly supported, 
for all t. 

Let us emphasize the two modifications with respect to the Definition in [22]: the first one, already done in [9], 
is that we use the spectral radius of the whole matrix d(i'), not only its negative part. The second one is that we 
deal with any a £ R, not only the case a = 0. 

The key property of this kind of solutions is that it allows the following weak-strong uniqueness result: if u 
is a smooth solution of the incompressible Euler equation (|70 p -(|7i p with uq as initial data then u = u. This can 
be easily obtained by observing that A{u) (respectively u — u) being a gradient (resp. divergence free) then the 

integral / A(u)(u — u)dx yamshcs. 

Theorem 6. Let be given a e M and uq G Ea- Let be given some smooth solutions [io^^f^)^, with a nice decay of 
/" as ^ goes to oo, of the equations (|67p -(|68 p - (|69p corresponding to some smooth compactly supported initial data 
[uJq, f^)s such that 



(wq)^ is bounded in L (R ) and p^ / f^dn is bounded in L'{W) (75) 



^This means d{v)ij := ^{djVi + diVj). 
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and such that, when e — ^ 0"*", 



e f lei Vo Odxd^ ^0, / - uol^dx ^ 0, 



(76) 




where Uq := K[ujq + Pq]. Then, up to an extraction, the sequence (u'^)e converges in L°°((0,T); Ea — w) to a 
dissipative solution of the incompressible Euler equation with initial condition uq . 

Note that as a simple corollary of Theorem [S] and of the weak-strong uniqueness result mentioned above we 
have in the particular case where uq is smooth with bounded vorticity that the whole sequence {u'^)e converges to 
the unique smooth solution of the incompressible Euler equation with initial condition uq. Let us refer here again 
to [HI and [23] for the existence and uniqueness of smooth solutions of the incompressible Euler equations in the 
plane with infinite energy. 

Proof. Notice first that the sequence (u'^)e(i) takes values in Ea, since 



Consider a smooth (in time/space) vector field v such as v{t) G E^ and curlw(t) is compactly supported, for all t. 
Let us denote 



The first thing to say about the modulated energy Hl{t) is that it is the sum of two nonncgative finite terms. In 
order to explain in what sense it can be understood as a modulation of the energy 




^"(0, x) + ^""(0, x)dx = a. 



2niit):^e \^-v{t,x)\^f%t,x,^)dxdC+ \u%t,x) - v{t,x)\^dx. 



(77) 





{p' + uj'){t, x){p' + uj'){t, y)G{x - v)dxdy, 



(78) 



one can reformulate it in the following way. Let us denote h curlu so that 



V = K[h] = V^(G*/i). 



(79) 



Then, using ^ and we get 



2Ul{t)^e[ \^-vit,x)\^rit,x,Odxd^- f {p' + - h){t,x){p' + u:' - h){t,y)G{x - y)dxdy, (80) 




so that one recovers the energy W{t) by setting w = in ([5(1)) . 

It is also not difficult to deduce from the assumptions ([TS]) and ([76]) that 




(81) 



Moreover the dynamic of Jif, is given by the following lemma: 



Lemma 9. We have 




v)'^{t, x)dx 



v{t,x))f'{t,x,Odxd£_ 



(82) 
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Recall the following classical notations, with the repeated index convention, 

1. ^ (X) C is the matrix {S.iCj)i,j and, given A a matrix, V : A is the vector {diAij)j. 

2. li A and B are two matrices, A : B is the scalar AijBij. 

Proof. To prove this lemma we proceed as in the Appendix of [3], with a few modifications. Let us observe that 
for any time, f^it), ui^{t) have a compact support. Expanding the square of the right hand side of (|80p we get 

2Hl{t) = eM2{r-){t)+e I p''\v\'^{t,x)dx - 2e f f{t)-v{t)dx 

- / {p' -h)(t,x){p'' +uj' ~h)(t,y)G{x-y)dxdy, 

hence 

2'Hl{t)^2W{t)=e I p''\v\^{t,x)dx-2e 1 f {t) ■ v{t)dx ^ ( h{t, x)h{t,y)G{x ~ y)dxdy 
+ 2 {p^ + u}^){t,x)h{t,y)G{x - y)dxdy. 

jR2xR2 

Therefore using that the energy T-L^{t) is independent of the time, we get 

27i 2/2 2/3 2/4 



2ni'{t)=e I dtp^v\'^{t,x)dx + ej p''dt\v\'^{t,x)dx -2e j dtf{t)-v{t)dx -2e j f(t) ■ dtv{t)dx 

2/5 



dth{t, x)h{t, y)G{x — y)dxdy 

2/6 2/7 



+ 2/ {p' +uj'){t,x)dth{t,y)G{x~y)dxdy + 2 dtip' + uj'){t, x)h{t,y)G{x - y)dxdy . 

jR2xR2 jR2xR2 

We are going to deal first with and then with I3. By integrating (|68|) with respect to ^ we get 

atp^ + div,(/) =0. (83) 

Adding this to (|57|) then yields 

+ w^) + div;r[/ + uj^u^] = 0, 

We therefore have 

/7 = / V(G * h) ■ [f + w'w'^lda; = / v ■ [f + uj'u']^dx, 



by using ((79|) . 

Let us also observe that multiplying ([68]) by ^ and then integrating with respect to ^ implies the following 
vectorial equation 

dtf + ^ ■ ^ ® ^f"d^] = - Z''"')^- (84) 

Using this we get also 



h = e I V 



JR2 |^Jr2 j 



dx + p'^v ■ u'^ dx — V ■ dx 

JR2 JR2 
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And since V 



(p*^ + cj^) = — which divergence gives back p'^ + 0;*^, we have 



B2 



/ u^^ i^u"^ ■.Vvdx+ u'-'- -Vu"^ ■ vdx. 

7r2 JR2 



where the second integral vanishes since 2u^ ■ Vm'^ = V 
get by symmetry 



and V is divergence free. We eventually 



/3 + /7 = -£ / d(w) : M ^ ® ^/de rfa; + / d{v) : u^^ (g) M-'^dx, 
that we decompose into: 

Qi 



l3 + l7=-e[ d{v):\[ i^~v)(E}{C-v)fdS,\dx+ [ d{v) : {v^ ~ u'^) ® {v^ - u'^)dx 



Q2 



+ -£ / yv: f ®vdx + I \/v : V f dx + I Vw : (g) da; + / Vw : (g) u 



+ £ pVu : u (g u H / V?j : ® dx 

■/r2 JR2 



Therefore we have 



ni'it) = h + i2 + h + h + h + Qi + - + Qi- 



Let us see a first cancellation: we have, with the repeated index convention, 



-e / {djv.i)vif^dx = - I \v\ divx{f )dx = -h, 



by using ((83 
Thus 



Hl'it) - /2 + /4+/5 + /6 + Q1 + Q2 + O4 



Q7 



= Qi + {h + I4 + Q2 + Qe) + ih + I6 + Q4 + Q5 + Qi)- 

Let us observe that Qi is equal to the sum of the first two terms of the right hand side of ([5^ . so that, in order to 
prove Lemma [UJ it is sufficient to prove that the third term is equal to I2 + I4 + Q2 + Qe ^-iid that the fourth term 
is equal to 1^ + 1^ + Qa + Q5 + Qi (here we have distinguished the terms with a factor e and the other ones). 
Regarding the third term, we have 

A{v)-{p'v-f)dx = \ [ p'dt\v\^dx^ 



dtv-fdx + J Vj{djVi)p^Vidx ~ j Vj{djVi)j^dx 

1 / p'^dt\v\'^dx — I dtv-fdx+ j p'\/v : v ® v dx — j yv:f®vdx, 

2 Jr2 JR2 Jr2 Jr2 



and hence, multiplying by £, we get: 



e / A{v) ■ {p^v - f)dx ^ I2 + h + Q2 + Q6 

iR2 
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We finally have (we use divxV = for the first equality) 

1 



A(v) ■ (v — u^)dx 



2 



df\v\'^dx 



dfV ■ u'^dx 



Vj{djVi)ufdx 



Wv : u"^ (E)v dx, 



so that it only remains to prove that 

- / Vw : (g) w = Q4 + Qs + Q7. 

Actually we are going to prove that 

V-y : ® u da; + = and Q5 + Qt = 0. 

2 

In order to do so, let us first note that we have the following relations for a, b and c three vector fields 



c • Va 



= c-Va Va:6( 



c • Va 



• b ^ Va-L : b^ 



Hence 



Vw : u'^ ^ V dx + Q4 = 



Vv 



U ® V + V (X) M 



W^""" (g) V — V ® u'^'^ 



dx 
dx. 



and since v = K[h] = V^(G* h), we get 

Vu : (E) V dx + Qj^ 



djd^{Gi<h) 



U: Vi — ViU, 



dx — 0, 



by symmetry. 

On the other hand 



Q5+Q7 



dx 
dx 



djdi{G -kh){ul — Vi)dj{G -kWjdx 
didj{G * — Vi)dj{G -k h)dx 



-V(G*/i) 



(u^ — u)da; = 0. 



□ 



With Lemma [9] in hands the proof of Theorem |6] then follows the proof of Theorem 6.1 in [9], with again a few 
modifications. For instance a crucial estimate in the strategy of [9] is that ||-\/ej^||Lf is bounded uniformly in 
£. Unlike Brenier, we cannot use here the energy of the system as the second term in ([75]) does not have a definite 
sign. Instead we are going to bootstrap an estimate involving both ||-\/£j'^||L~(Li) and ||H^||loo. 

By Cauchy-Schwarz inequality we get 



Vi / \f\{t,x)dx ^ V~e I \^-vit,x)\rit,x,Odi 

^ {{ni{t))^ + \\vU^)u\\L^- 



dx 



\v(t, x)\p'(t, x) dx 



(85) 
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Above we used that (jSSll implies that ||/3'^(t)jlLi = UpoIU^ ^o^' ^■ 
Now, using lemma [SJ we deduce that 

ni\t)^2\\d{v)u^ni{t) + s\\A{v)u^\\vU^y^^^^ 

that is, 

K'it) ^P,,,[\\v\\wi.^ +\\v\\^^^^.,2^)[l+HtXt)], 

for some polynomial function Pin. Here we used the assumption ([75]) to bound ||poIIli- 

A Gronwall lemma gives then that for a fixed smooth vector field v (such as v{t) S Ea for all t), 



HI is bounded inij", 
which implies, according to (|77p and (|55|) . that 

{Vsf)e is bounded in 2.^(^1^) 
{u" - v)e is bounded in Lf{Ll). 

Let us now look for a bound of p'u' . We start by writing 

K[p^+ui^-h] 

p^u^ = {p^ +uj^ - h) 'i^F^^ +{h-uj^){u'' -v)+ p'v. 



(86) 

(87) 
(88) 

(89) 



From ([57]) wc infer that ||a;^(t)||^2 = ||a;o||L2 for any t. Using the assumption ((75|) we deduce that {u!'^)e is bounded 
in L'^{L1), and we recall that (p'^)e is bounded in Lf[Ll.). Since /i,?; g L^, using (|88| we see that the two last 
factors of the right hand side are bounded in L'f°{Ll.). 

For the first one we will use, as in [9], the following lemma. 



Lemma 10. For any g smooth such as / g{x)dx = 0, we have 

JR2 



\9K\ 



(90) 



Let us provide a proof for sake of completeness. 
Proof. Let if he a, smooth compactly supported vector field over R^. We have 

ip{x) ■ g{x)K[g\-^{x)dx = 







V 




dx 


/R2 


G^g 


a, 


G^g 


dx 



ipid-j 



G*g 



d-jdi 



G^g 



dx 



G^g 



G*g 



div((^)(x)- 



VG*.g 



G-kg ■ W{ip) Girg 



dx 

dx. 



and hence 



/ 



'Pix) ■ gK[g]{x)dx 



^ -M,,n..^\\K[g]\\L^, 



what yields the result. 



□ 
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The lemma above thus gives that the first term in the right hand side of (|89p is bounded in L'^(W~^'^), so 
that p^u'^ is eventually bounded in for some m large enough. 

Now we are going to use this to deduce that {p^)^ is relatively compact in Wc first combine 

and (HH) to get 



p'u'. 



Because of fSS]) , the first term of the right hand side is bounded in '£^{yV,j, ^'^) that we already injected in some 
in which the second term is also bounded. Since H-"'{B{0,N)) ^ ff-™-i(B(0, iV)) is compact for 
all N € N, Ascoli's theorem and a diagonal extraction allow us to conclude that (p"^ — sdivx{j^'^))e is (up to an 
extraction) converging in '^^{^'{R'^)). But {y/ef)^ is bounded in %^{Ll) ^ '^i°(^'(R^)) and hence we get the 
convergence of (p'^)e in '^°(^'(R2)) and we denote by p its limit. 

The previous convergence holds (up to an extraction) as well for {uj^)e i-^, the reasoning being simplified 
because (a;'^)^ is bounded in every L'^{LP). Wc hence get a subsequence of {u^)e ~ {K[p^ +uj^])^ that converges in 
^t°(^'(R2)) to some u. But {u" - v)e is bounded in '^t°(i^) so that the previous convergence implies (u^ — v)^ — > 
(u — v) in "t^^ (L^. — w) , with furthermore the pointwise bound 

\\u{t)-v{t)\\L2 i^\jm\\u,{t)~vit)\\L2. 

£->0 

(HJ)e is up to an extraction converging in — w-k to some function and the previous bound gives 



almost for all t G [0,r]. Now lemma [5] gives, 

nl'it) ^2\\div)\\it)nlit) + [ Aiv){t,x) ■ iv-u''){t,x)dx 
hence 

•H,^(0 s;H^(0)exp( / 2\\d{v)\\{s) ds\ + f cxp ( / 2\\d{v)\\{a)da} I A{v){s,x) ■ {v - u'){s,x)dxds 



(91) 



1 



\v\\w^. 



4r\\LriLl)+4Pl\L^{Ll) 



1 



4f\\L^{Ll) +£\\p^\\l^{li) 



Using previous bounds, we see that the product of the second line goes to with e. We now pass to the limit in 
the previous inequality, using (|8ip . to obtain 

■Ht,(t) =^ ||iio-w|||2Cxp|^ 2||d(w)||(.s)ds| + J cxp 1^ 2||d(u)||(cr)dcr| y A{v){s,x) ■ {v - u){s,x)dxds, 

first in the distributional sense (that is against any positive test function) and then almost everywhere since both 
sides of the previous inequality are L°° functions. 

Then (|9ip allows to obtain the dissipation formulation introduced in Definition [6l at least for any smooth (in 
time/space) vector field v such as v{t) S and compactly supported curlz;(t), for all t, which is sufficient in view 
of Remark HI □ 
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